Measurement of the tt Production Cross Section in pp Collisions at ^/s = 1.8 TeV 
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(February 7, 2008) 

We present a measurement of the ti production cross section in pp collisions at ^/s = 1.8 TeV 
using an integrated luminosity of 109 pb~^ collected with the Collider Detector at Fermilab. The 

measurement uses ti decays into final states which contain one or two high transverse momentum 
leptons and multiple jets, and final states which contain only jets. Using acceptances appropriate 
for a top quark mass of 175 GeV/c^, we find cr^j = 7.61J J pb. 

14.65 Ha, 13.85 Ni, 13.85 Qk 
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The measurement of the pp ^ tt X production cross section presents a test of both the production and decay 
mechanisms of the standard model. Recent calculations based on Quantum Chromodynamics (QCD) have led to 
predictions for the cross section with a theoretical uncertainty of less than 15%. A measurement that is significantly 
different from the predicted value can signal either non-standard model production, for instance the decay of a 
heavy resonance into tt pairs, or a non-standard model decay mechanism such as the decay into supersymmetric 
particles In the latter case it is of particular interest to measure the cross section into different final states, 
because an unexpected decay mode of the top quark will modify the expected branching fractions. The tt production 
cross section has been measured before by both the Collider Detector at Fermilab (CDF) and DO collaborations 

The standard model predicts that the top quark will decay nearly 100% of the time to Wb. The W boson can then 
decay to either a pair of quarks, or a lepton neutrino (iv) pair. We categorize the decays of ti pairs by the decays 
of the two W bosons as either lepton-|-jets, dilepton, or all-hadronic. The dilepton and all-hadronic analyses are 
described elsewhere . We now have nearly twice as much data as reported in |Q . With improved measurements of 
acceptances and backgrounds, and by combining all the decay modes, we measure the cross section with better than 
twice the precision of our previous measurement. 

The data presented here represent the entire data set accumulated between 1992 and 1995 with the CDF detector, 
and corresponds to an integrated luminosity of 109±7 pb~^ (19 pb~^ from the 1992-93 run and 90 pb^^ from the 
1994-95 run) §]. 

The CDF detector consists of a magnetic spectrometer surrounded by calorimeters and muon chambers. A four- 
layer silicon vertex detector (SVX), located immediately outside the beam pipe, provides precise track reconstruction 
in the plane transverse to the beam and is used to identify secondary vertices from b and c quark decays. A detailed 
description of the detector can be found elsewhere . 

The electron, muon, and multi-jet events used in this analysis were selected by a three-level trigger. Lepton samples 
were acquired with inclusive electron and muon triggers requiring Pr(lepton)> 18 GeV/c. A missing transverse 
energy trigger was also used in order to recover events lost due to small inefficiencies in the inclusive lepton 

triggers. 

Decays of tt pairs into lepton-f jets are characterized by a single high-Pr lepton, missing transverse energy from the 
W ^ Hv decay, plus four jets, two from the hadronically decaying W boson and two from the b quarks from the top 
decays. Jets are defined using a cone algorithm with Ai? = \/ At/)^ + /S.rf—OA, where 77 is the pseudo-rapidity. Jets 
are counted in this analysis if | 77 |< 2.0. The number of observed jets may decrease due to detector effects or jet 
overlap, or increase as a result of multiple interactions or the presence of gluon radiation. In the lepton-f jets channel, 
events with three or more jets with measured Ft >15 GeV define the tt signal region. 

The data sample for the lepton-|-jets analysis is a subset of a sample of high-Pj. inclusive lepton events that contain 
either an isolated electron with > 20 GeV or an isolated muon with P,, > 20 GeV/c in the central region 
(I 77 |<1.0). Events that contain a second same flavor lepton of opposite charge are removed as Z boson candidates 
if the reconstructed ee or jjLjji invariant mass is between 75 and 105 GeV/c^. If a candidate high-ET photon |]l^ is 
present, the three-body mass is used to remove radiative Z candidates. An inclusive W boson sample is selected from 
the inclusive lepton sample by requiring > 20 GeV and that the lepton be isolated from any jet activity. For the 
latter we define isolation, /, as the transverse energy in a cone of Ai? = 0.4 centered on the lepton, but excluding the 
lepton energy, divided by the Ft (Ft) of the electron (muon), and require / < 0.1 . Furthermore, the event must not 
be accepted as a dilepton candidate 

In order to separate ti events in the lepton+jets channel from the large W^-l-jets background, we require that one of 
the jets be identified as a 6 jet candidate. Identification of b jets is done either by reconstructing secondary vertices 
from 6-quark decay using the SVX (SVX tagging), or by identifying an additional lepton from a semileptonic b decay 
(SLT tagging). The SVX and SLT tagging algorithms are described in Ref. 

The efficiency for tagging a b quark from a top decay is determined from tt Monte Carlo data together with a 
detailed simulation of the detector, which includes the effects of local track density on the track finding efficiency. 
The systematic uncertainty due to the tracking efficiency modeling is determined by comparing data and Monte Carlo 
tracking efficiencies and multiplicity distributions as a function of jet E^ in inclusive electron and muon samples, 
which are enriched in b decays. The efficiency for tagging at least one b quark in a ti event with > 3 jets is found to 
be (39±3)%. Of this 39%, a factor of 67% comes from the fiducial acceptance of the SVX. 

The SLT algorithm identifies both muons and electrons with P,, > 2.0 GeV/c to | 77 |= 1.0. The efficiency of this 
algorithm, as a function of lepton Pr, is measured with photon conversion and J/t/j /i/i data, and applied to Monte 
Carlo ti events. The probability of finding an additional e or from a b quark decay in a ti event with > 3 jets is (18 
± 2)%. 

In the tt signal region of W+ > 3 jets, there are 34 SVX-tagged events containing a total of 42 SVX tags, and 
40 SLT-tagged events containing a total of 44 SLT tags. Of these, 11 events are tagged by both the SVX and SLT 
algorithms. 
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The acceptance for identifying tt events in the lepton+jets mode is calculated from a combination of data and 
PYTHIA jl^ and HERWIG @ tt Monte Carlo samples. We use a top mass of 175 GeV/c^ |jl^ when evaluating the 
acceptance. The total acceptance, including the branching fraction, is calculated as the product of the kinematic 
(including lepton identification) and geometric acceptances, the trigger efficiency and the tagging efficiency. We 
measure these efficiencies as described in Rcf. |3|, and average over the two running periods. For the lepton+jets 
analysis, the product of the geometric and kinematic acceptance is (10.4±1.0)%, and the trigger efficiency is (90±7)%. 
These factors are common between the SVX and SLT analyses. Combining with the respective tagging efficiencies 
gives a total SVX acceptance of (3.7±0.5)% and a total SLT acceptance of (1.7±0.3)%. 

The systematic uncertainties on the geometric and kinematic acceptances come from the following: the jet energy 
scale (±5%), modeling of initial state gluon radiation (±2%), final state gluon radiation (±5%), Monte Carlo depen- 
dence and modeling (±5%), detector resolution effects (±2%) and instantaneous luminosity dependence (±1%). The 
uncertainties on the tagging and trigger efficiencies are dominated by the level of agreement between data and the 
Monte Carlo predictions. 

The most important source of background in the SVX-tagged lepton+jets channel is inclusive W production in 
association with jets containing b or c quarks, eg. pp — > Wg (g — *■ 66). In addition, there are contributions to the 
background from mistags (i.e. tags in jets which contain no true displaced vertices), and small contributions from the 
following processes: non-VF(e.g. direct 66 production), single top production, WW, WZ, and Drell-Yan. 

To calculate the background from VF+heavy flavor events, we use the herwig and vecbos jlj] Monte Carlo 
programs to predict, as a function of jet multiplicity, the fraction of ly+jet events which are Wbb, Wcc and Wc. 
These fractions, and a tagging efficiency for each type of event, are applied to the number of VF+jet events seen in 
the data to give an expected background from these sources for each jet multiplicity. The details of this method can 
be found in Ref. [§. 

The background from events in the sample that do not contain real W bosons (non-M^) is calculated from the data 
by measuring the number of tags as a function of lepton isolation, /, and Ei^- The tagging rate in the low EJt, high-/ 
region, where there are essentially no real W events, is used to predict the contamination in the W signal region of 
high EJt, low /. 

To calculate the background from mistags we assume that the distribution of reconstructed transverse decay 
length, Lxj,, from this source is symmetric about zero. Secondary vertices with negative h^y (i.e. those which 
reconstruct to the opposite side of the primary from the jet direction) come primarily from reconstruction errors in 
light quark jets. We parametrize the negative L^j, distribution measured in generic jet data as a function of jet Et, rj, 
and the number of SVX tracks in the jet. This parametrization is applied to the W^+jets data to predict the number 
of mistags observed. 

The single top background is determined by measuring the acceptance for W* and W^-gluon production using the 
PYTHIA and HERWIG Montc Carlo programs, and using the latest theoretical cross sections jl^. The remaining, 
relatively small, backgrounds {WW, WZ, Z tt) are derived from Monte Carlo predictions. 

The individual components of the background and their totals are shown in Tabic |. In addition to the signal region 
of 3 or more jets, we show the predicted number of tags in events with 1 and 2 jets as a check of our calculation. 
An iterative correction, to account for the tt content of the VK+jet events, is applied to those backgrounds that are 
calculated as a fraction of the observed number of these events g . The corrected background in the signal region is 
9.2+1.5 tagged events. We observe 34 tagged events, resulting in a cross section of 6.2't_i'\ pb. 

The background to SLT-tagged events is dominated by W events with hadrons misidentified as leptons (including 
decays in flight), with electrons from unidentified photon conversions, or with real heavy flavor jets {Wbb, Wcc). 
These backgrounds are calculated by measuring the fraction of tags per track in a generic jet sample as a function of 
the track Ep. These fractions are applied to tracks in the M^+jet events to estimate the background from the above 
sources. Smaller backgrounds are, in order of importance in the signal region, WW and WZ, non-M^, Z tt, single 
top, Wc, and Drell-Yan production. The results of the background calculation and the number of tags observed in 
the data are shown in Table |. In the signal region of W+ >3 jets, the background prediction is 22.6 + 2.8 tagged 
events. We observe 40 SLT tagged events, resulting in a cross section of 9.2 Ig g pb. 





W+1 Jet 


W+2 Jets 


W+3 Jets 


W+ >4 Jets 


Events before tagging 


10 716 


1663 


254 


70 


SVX tagged events 


70 


45 


18 


16 


Wbb + Wcc 


19.3 ±6.7 


9.7 ±2.4 


2.3 ±0.6 


0.85 ±0.24 


Non-W 


7.7 ±3.0 


4.0 ±1.5 


1.4 ±0.5 


0.77 ±0.33 


Mistags 


20.9 ±6.3 


7.2 ±2.1 


1.7 ±0.5 


0.63 ±0.22 


Single top 


1.3 ±0.4 


2.8 ±0.7 


1.0 ±0.4 


0.29 ±0.14 
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SLT tagged events 
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78 


25 


15 


Mista.gs~^W bb~\-Wcc 


235 ±21 


66.6 ±6.1 


15.1 ±1.4 


6.8 ±0.7 


Single top 


0.9 ±0.3 


1.5 ±0.5 


0.6 ±0.3 


0.2 ±0.1 


Other 


33.1 ±10.6 


9.6 ±3.0 


1.2 ±1.4 


0.6 ±0.6 


Uncorrected bkgnd total 


269±23 


77.7±6.6 


16.9±2.0 


7.6±0.9 


Corrected bkgnd total 


269±23 


77.7±6.6 


15.9±2.0 


6.7±0.8 


SLT tagged tt expected 


0.8±0.2 


3.8±1.2 


6.6±1.7 


7.7 ±2.1 



TABLE I. Summary of event yields from the lepton±jets analyses. The expected tt contribution is calculated using the 
measured combined cross section from this paper. 

Our best measurement of the it cross section comes from combining the results of the lepton+jets analyses with 
the dilepton and all-hadronic analyses . The results of the individual analyses are summarized in Table The 
dilepton analysis finds 9 candidate events, with an expected background of 2.4±0.5. The all-hadronic analysis has two 
parallel paths, one which requires a single SVX tagged jet plus kinematic cuts to isolate tt^ and a second which requires 
two SVX tagged jets, but no additional kinematic cuts. The single tag analysis identifies 187 candidate events with an 
expected background of 142±12 events, while the double tag analysis identifies 157 candidates and predicts 120±18 
background events. There are 34 candidate events in common between the two analyses. The dilepton, lepton±jets, 
and all-hadronic data samples are exclusive sets. 

We calculate the tt production cross section from the combined results of the dilepton and lepton+jets channels using 
the same maximum likelihood technique described in ^ . The all-hadronic result is added by including the multivariate 
Gaussian term described in Ref. 0. The likelihood properly accounts for correlated systematic uncertainties, such 
as the uncertainty on the integrated luminosity, and the uncertainty on the lepton±jets geometric and kinematic 
acceptance, which is common to both the SVX and SLT analyses. 





Lepton+Jets 


Dilepton 


All-Hadronic 


Tag 


SVX 


SLT 


not req. 


SVX 


2 SVX 


etag 

^trigger 


0.39±0.03 1 0.18±0.02 
0.104±0.010 
0.90 ±0.07 


0.0076±0.0008 
0.98±0.01 


0.42±0.04 
0.106±0.021 



0.11±0.02 

0.263±0.045 
qqQ+0.002 1 


Etotal 


0.037±0.005 


0.017±0.003 


0.0074±0.0008 


0.044±0.010 


0.030±0.010 


Obs. Events 
Background 


34 
9.2 ±1.5 


40 

22.6 ±2.8 


9 

2.4 ±0.5 


187 
142 ±12 


157 
120 ±18 


o-ti (pb) 


6.2l^:^ 


9 2 +^'^ 

^■^ -3.6 


8-211:1 


9 fi+'* " 

^•"-3.6 


n 5+''''' 



TABLE II. Summary of acceptance factors and measured cross sections for each analysis channel. The acceptances are 
calculated for a top quark mass of 175 GeV/c^. 
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In Figure |i| we show the results of the cross section calculation for each tt decay channel, as well as the combined 
measurement. The combined cross section for Mtop = 175 GeV/c^ is T.et^J pb , where the quoted uncertainty 
includes both statistical (±1.2 pb) and systematic effects. Due to the mass dependence of the acceptances, the 
calculated cross section changes by ±10% for a ±15 GeV/c^ change in top mass from 175 GeV/c^. Theoretical 
calculations [Q| range from 4.75 pb to 5.5 pb for Mtop = 175 GeV/c^. From the ratio of the measured cross sections in 
the dilepton, lepton+jets and all-hadronic channels we can calculate the branching fraction for a top quark decay to a 
final state electron or muon. Assuming lepton universality and W decay acceptance, the apparent branching fraction 
to an electron or muon is 0.188+/-0.048, consistent with the standard model expectation of |. Specifics of possible 
non-standard model top decays have not yet been considered. 
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FIG. 1. Measured ti production cross section at Mtop=175 GeV/c'^ for each of the decay channels and for the combined 
measurement. The lepton+jets cross section is calculated from the SVX and SLT analyses described in the text. The line is 
drawn through the central value of the combined measurement. The theory point shows the spread in the central values of the 
3 most recent predictions 
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